Lattice Index Theorem and Fractional Topological Charge 
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We study topological properties of classical spherical center vortices with the low-lying eigenmodes 
of the Dirac operator in the fundamental and adjoint representations using both the overlap and 
asqtad staggered fermion formulations. In particular we address the puzzle raised in a previous 
work of our group [Phys. Rev. D 77, 14515 (2008)], where we found a violation of the lattice index 
theorem with the overlap Dirac operator in the fundamental representation even for "admissible" 
gauge fields. Here we confirm the discrepancy between the topological charge and the index of the 
Dirac operator also for asqtad staggered fermions and the adjoint representation. Numerically, the 
discrepancy equals the sum of the winding numbers of the spheres when they are regarded as maps 
R 3 U {00} — > SU(2). Furthermore we find some evidence for fractional topological charge during 
cooling the spherical center vortex on a 40 x 2 lattice. The object with topological charge Q = 1/2 
we identify as a Dirac monopole with a gauge field fading away at large distances. Therefore even 
for periodic boundary conditions it does not need an antimonopole. 
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I. INTRODUCTION 

Topological models have been introduced to explain 
quark confinement in non-abelian gauge theories. There 
are good reasons to believe that a force strong enough 
to conhne quarks must also break chiral symmetry spon- 
taneously P|. The center vortex model 0t3| seems to 
be the most promising candidate to explain both phe- 
nomena. Center vortices, quantized magnetic flux lines, 
compress the gluonic flux into tubes and cause a lin- 
early rising potential at large separations. Numerical 
simulations have indicated that vortices could also ac- 
count for phenomena related to chiral symmetry, such as 
causing topological charge fluctuations and spontaneous 
chiral symmetry breaking (SCSB) fa-tllT]. These non- 
perturbative features of the QCD vacuum are intimately 
linked to the properties of the low-lying spectrum of the 
Dirac operator. In Ref. [12| we found strong correla- 
tions between the center vortices and the density distri- 
bution of low-lying Dirac eigenmodes. We showed that 
even center-projected lattice configurations give rise to a 
spectral density of near-zero modes, which by the Banks- 
Casher relation [13] is proportional to the chiral conden- 
sate, the order parameter for SCSB. Another relation 
between fermions and to polo gy is given by the Atiyah- 
Singer index theorem (l4H17| . It states that the topolog- 
ical charge of a gauge field configuration equals the index 
of the Dirac operator in this gauge field background. In 
Ref. 18[ we investigated the lattice index theorem and 



the localization of the zeromodes for thick classical center 
vortices. For nonorientable spherical vortices, the index 
of the overlap Dirac operator turned out to differ from 
the topological charge. This may be related to the fact 
that even in Landau gauge some links of this configura- 
tion are close to the nontrivial center elements. 

In this paper, we work with thick spherical vortices 
in SU(2) lattice gauge theory and extend our analysis 
of the topological charge and the lattice index theorem. 



In Sec. [H] we re-introduce the construction of a spheri- 
cal vortex and report on our calculations with the over- 
lap Dirac operator. After analyzing the response of (im- 
proved) staggered fermions to the topological background 
of the gauge field we investigate in Sec. IIIII the localiza- 
tion of asqtad staggered zeromodes with respect to the 
position of the thick vortices and find the same discrep- 
ancy for the topological charge determined by different 
methods. We then examine the index theorem for ad- 
joint fermions in Sec. IIVI and discuss the role of adjoint 
fermions with respect to fractional topological charge in 
Sec. IV] We find an interesting configuration giving topo- 
logical charge one half which we analyze in Sec. IVII in 
more detail. We end with our conclusions. 



II. SPHERICAL CENTER VORTICES AND 
DIRAC ZERO MODES 

The nonorientable spherical vortex of radius R and 
thickness A is constructed with the following links: 
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where r is the spatial part of x v and the profile function 
a is either one of a+, a-, which are defined by 
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This means that all links are equal to 1 except for the 
t-links in a single time-slice at fixed t = 1. The phase 
changes from to it from inside to outside (or vice versa) . 
The grap h of a_ (r) for a 40 3 x 7V t -lattice is plotted in Fig. 
2 in [18j |. giving a hedgehog- like configuration, without 
any intersections and hence no topological charge. Since 
only links in the time direction are different from 1, the 
topological charge determined from any lattice version of 
FF vanishes for this spherical vortex configuration. We 
shall refer to such determinations as topological charge 
in the rest of this paper. 

However, we find a discrepancy to the index of the over- 
lap Dirac operator [lj, LLil, [20] . According to the Atiyah- 
Singer index theorem the topological charge is related to 
the index by ind D[A] = rt_ — n + = Q where n_ and 
n + are the number of left- and right-handed zeromodcs 
of the Dirac operator |14| - KL6| . For the spherical vortex 
configuration, Eq. (J2.2I) . we obtain a nonzero index. 

The lattice version of the index theorem is only valid 
as long as the gauge field is smooth enough and satisfies 
a so-called "admissibility" condition. This condition as- 
sures that HJ has no zero eigenvalues so that the sign 
function in the definition of the overlap Dirac operator 
is well defined. It requires that the plaquette values U^ v 
are bounded close to 1, the value for very smooth gauge 
fields. A sufficient, but not necessary bound for the "ad- 
missibility" of the gauge field is tr(l-C/ Aliy ) < 0.03 [H|-|3 

On a 40 3 x iV t -lattice the traces of the plaquettes for the 
spherical vortex, Eq. (12.21) . deviate only by a maximum 
of 1.5% from trivial plaquettes, satisfying the "admissi- 
bility" condition. Nevertheless, the index of the overlap 
operator is nonzero, ind D = =pl, for a±. 

To try to understand the discrepancy we apply stan- 
dard cooling to the spherical vortex configuration. For 
many cooling steps, the index of the overlap Dirac oper- 
ator does not change, but the topological charge quickly 
rises close to =pl for a± while the action S reaches a 
(nonzero) plateau, as shown in Fig. Q] for a 40 3 x 4 lat- 
tice. So, the index of the overlap Dirac operator agrees 
with the topological charge after some cooling, but not 
on the original vortex configuration. 



III. TOPOLOGY AND STAGGERED 
FERMIONS 

Staggered fermions don't have exact zeromodes, but 
a separation between "would-be" zeromodes and nonchi- 
ral modes is observed for improved staggered quark ac- 
tions [24]]. These results shall be verified first, by com- 
paring standard and asqtad improved staggered fermions 
and by u sing the HYP (hypercubic blocking) smearing al- 
gorithm [251 ] . The results are shown in Fig. [5] presenting 
the first twenty eigenmodes of thirty configurations, gen- 
erated by lattice Monte Carlo simulation of the tadpole 
improved Liischer-Weisz pure-gauge action at coupling 
Plw = 3.7 (lattice spacing a « 0.09 fm). We chose a 
fairly small lattice spacing to keep lattice effects as small 
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1. Cooling of a spherical vortex on a 40 3 x 4 lattice. The 
topological charge rises from zero to close to one for a — a_ 
(right scale) while the action S (in units of the one-instanton 
action Sinst) reaches a plateau (left scale). 



as practicable. The eigenvalues are plotted against the 
chirality (pseudoscalar density) of the modes, which for 
staggered fermions is determined by ('JTs'F), where T$ 
corresponds to a displacement along the diagonal of a hy- 
percube. To ensure gauge invariance the product includes 
gauge field multiplications along all shortest paths con- 
necting opposite corners of the hypercube. Comparing 
standard (red + symbols) and asqtad improved (green 
x symbols) staggered fermions, a slight improvement is 
observed for the latter: "would-be" zeromodes show a 
higher chirality (up to 0.5) and eigenvalues are "closer" 
to zero. But for these still rather rough configurations 
it seems hard to really identify the would-be zeromodes 
reliably. To make the configurations smoother, while re- 
taining topological fluctuations, we applied five steps of 
HYP-smearing. Then the chirality of the would-be zero- 
modes gets well defined and the eigenvalues get close to 
zero. On these smoothed configurations there is no big 
difference between standard (blue * symbols) and asqtad 
improved (magenta □ symbols) fermions anymore. 



standard staggered fermions 
asqtad improved staggered fermions 

after 5 steps of HYP-smearing: 

standard staggered fermions 
asqtad improved staggered fermions 




0.15 
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FIG. 2. The first twenty eigenmodes of thirty Monte Carlo 
configurations on a 20 4 lattice, for standard and asqtad stag- 
gered fermions on the original and HYP-smeared gauge fields, 
with antiperiodic boundary conditions. 




FIG. 3. Scalar densities of asqtad staggered zeromode of pos- 
itive chirality for the positive (a+, left) and negative chirality 
for the negative (a_ , right) spherical vortex in a time slice of 
a 40 3 x 8 lattice, with antiperiodic boundary conditions. 



These results promise that one should be able to iden- 
tify zeromodes reliably on smooth configurations, such 
as the above examples of spherical vortices. The asq- 
tad staggered modes were calculated for a positive and 
a negative non-orientable spherical vortex made up of 
time links in a single time slice of a 40 3 x 8 lattice. 
Due to charge conjugation, staggered eigenmodes are 
doubly degenerated but according to the index theorem 
(ind D [A] = n_ — n + = 2Q) we get exactly the same 
topological charge as the overlap Dirac operator. The 
scalar densities of the would-be zeromodes for both types 
of spherical vortices show a distribution (Fig. [3]) that 
looks very similar to their overlap counterparts (see [18j ) . 



IV. INDEX THEOREM AND ADJOINT 
FERMIONS 

It is also interesting to study the correlation of the 
number of zeromodes and the topological charge for 
fermions in the adjoint representation, as discussed at 
recent lattice conferences (see e.g., j26|). The index of 
the massless Dirac operator in the adjoint representation 
of the SU(N) gauge group in a background field of topo- 
logical charge Q is equal to 2NQ [27j . We test the lattice 
index theorem for overlap and asqtad staggered fermions 
in the adjoint color representation on the spherical vortex 
configurations, Eq. (2.1). Since the fermion is in the real 
representation, the spectrum of the adjoint Dirac opera- 
tor is doubly degenerate. Therefore, the index can only 
be even valued. It will be interesting whether we find all 
possible even values or only multiples of four, resulting in 
a gauge field background made up of classical instantons. 

In Table U the number of positive and negative over- 
lap and staggered zeromodes in the fundamental and ad- 
joint representation for periodic and antiperiodic bound- 
ary conditions on positive and negative spherical vortex 
configurations on different lattice sizes are summarized. 
Concerning this table we like to mention two issues. All 
investigated fermion representations equally violate the 
lattice index theorem for these special configurations: the 
index is nonzero while the topological charge vanishes. 
Further, in some cases, indicated in bold, the adjoint in- 
dex is not a multiple of four, as will be discussed in the 



next sections. 

To study the above mentioned violation, we consider 
the homotopy of the above spherical vortex configura- 
tions, Eq. (2.1). The t-links of these spherical vortices 
fix the holonomy of the time-like loops, defining a map 
Ut(x,t = 1) from the xyz-hyperplane at t = 1 to SU(2). 
Because of the periodic boundary conditions, the time- 
slice has the topology of a 3-torus. But, actually, we can 
identify all points in the exterior of the 3 dimensional 
sphere since the links there are trivial. Thus the topol- 
ogy of the time-slice is R 3 U {oo} which is homeomorphic 
to S 3 . A map S 3 — > SU(2) is characterised by a winding 
number 



N=- 



1 



24tt 2 



d 3 x e. 



jk tx[(U''d i U)(U^d j U){U''d k U)}, 

(4.1) 



resulting in J\f = — 1 for positive and JV = + 1 for negative 
spherical vortices. With this assignment the index of the 
Dirac operator and the topological charge after cooling 
coincide with this winding number Af. Obviously such 
windings, given by the holonomy of the time-like loops of 
the spherical vortex, influence the index theorem [28l |29j . 
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2+0- 
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TABLE I. Number of positive and negative overlap and stag- 
gered zeromodes in the fundamental and adjoint represen- 
tation for periodic and antiperiodic boundary conditions on 
positive and negative spherical vortices on different lattice 
sizes. Numbers which indicate half-integer topological charge 
are printed in bold. 



V. ADJOINT ZEROMODES AND 
FRACTIONAL TOPOLOGICAL CHARGE 

Classical instantons carry an integer topological 
charge. Thus, in case of a fermion in the fundamental 
representation of SU(2) there is exactly one zeromode 
for a one-instanton configuration. Now, if the actual 
constituents of the QCD vacuum had topological charge 
Q = 1/2, no zeromode would be produced. However, for 
adjoint fermions configurations with topological charge 
Q = 1/2 are able to create a zeromode. Edwards et al. 
presented in [27| some evidence for fractional topologi- 
cal charge on the lattice. Garcia-Perez et al. (26|, how- 
ever, associated this to lattice artefacts, i.e., topological 
objects of size of the order of the lattice spacing. For 
comparison we repeat the calculations with Monte Carlo 
configurations on 12 4 lattices for overlap and 16 4 lattices 
for asqtad staggered fermions. In the first case we find 
for five out of twenty configurations a number of overlap 
zeromodes which is not a multiple of four, i.e., 2, 6, 6, 
10, 10. For adjoint asqtad staggered fermions multiples 
of eight are expected, since staggered fermions have a 
doubly degenerate eigenvalue spectrum. Nevertheless we 
even find a couple of configurations with one or two "ze- 
romodes" , but of course it is not clear whether all these 
low-lying modes are would-be zeromodes or whether all 
would-be zeromodes were identified. 

However, the measurements with adjoint fermions on 
the classical spherical vortices also show fractional topo- 
logical charge for small lattice sizes, see Table UJ Very 
interesting is the case of the cooled 40 3 x 2 lattice con- 
figuration, which will be discussed next. 



VI. COOLING OF VORTEX CONFIGURATIONS 

During cooling of a spherical vortex on a 40 3 x 2 lat- 
tice we find evidence for fractional topological charge as 
shown in Fig. @] It exhibits a second "plateau" for the 
topological charge with Q « 1/2 between cooling steps 
100 and 130. We analyzed this configuration with adjoint 
eigenmodes, which indeed measure fractional topological 
charge, as seen in Table HI The spherical vortex con- 
tracts during cooling as shown in Fig. [5j After 78 cool- 
ing steps the vortex structure vanishes. In the maximal 
abelian gauge one can identify a monopole-antimonopole 
ring after projection which again contracts during cool- 
ing and disappears after 91 cooling steps, see Fig. [B] In 
the interesting region with topological charge Q = 1/2 
(at cooling step 120), the action and topological charge 
densities concentrate in the center of the spacial volumes 
and the gauge field at the lattice boundary is trivial. Lan- 
dau gauge yields a very symmetric configuration, a static, 
singular monopole without an antimonopole. In Lan- 
dau gauge both time slices are exactly the same and all 
time-like plaquettes are trivial. The central cube in ev- 
ery time-slice is a non-abelian representation [30j, |3l| of 
a Dirac monopole, see also Fig. 9. Its six plaquettes 
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FIG. 4. Cooling of a spherical vortex with a — a~ on a 40 3 x 2 
lattice. While the action S (in units of the one-instanton ac- 
tion Sinst) decreases slowly (left scale), the topological charge 
first rises from zero to close to one (right scale), then decreases 
to an intermediate plateau of Q « 0.4 before it vanishes. 





cooling steps 20 cooling steps 40 cooling steps 




60 cooling steps 70 cooling steps 78 cooling steps 

FIG. 5. Cooling of a spherical vortex on a 40 3 x 2 lattices. 
The vortex shrinks and vanishes after 78 cooling steps. 



correspond to rotations by 2-7r/3 in the fundamental rep- 
resentation of SU(2). In the color frame of the corner 
with the smallest coordinate values the plaquette color 
vectors, parallel-transported according to Fig. 2 in [32J, 
point in the (— 1, — 1, — Indirection. Hence, the six pla- 
quettes sum up to a total rotation of 47r in agreement with 
the non-abelian lattice Bianchi identity, which states that 
the product of plaquette rotations in appropriate order 
results in the unit matrix. Nevertheless the magnetic flux 
out of the cube is nonvanishing [32] ■ 

In the U(l) representation in color direction (1,1,1) 
the central cube represents a Dirac monopole with pla- 
quette values summing up to 27r. The link variables of 
the central cube correspond to an SU(2) rotation of a 
unit vector from one corner to the other, i.e., cos(oj) = 
1/3 = (— 1, — 1, — 1)(1, — 1, — 1)/3. Due to the non-abelian 
nature of the SU{2) gauge field for increasing distance 
from the center the field strength approaches zero. No 
antimonopole is needed to compensate the monopole. 
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FIG. 6. Cooling of a spherical vortex with a — q_ on a 40 3 x 2 
lattice. In maximal abelian gauge a monopole-antimonopole 
circle shrinks and vanishes at 92 cooling steps. 
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FIG. 8. Polyakov-line at spacial site (20,20,20) for the nega- 
tive spherical vortex during cooling shows two transitions as 
function of the cooling step. 



> 

o 

a 
>. 





' ; •' 


sill 


M M-' 








08 




" . . " n » x . 






t »• 


» » " - .■%"* + 




0.6 


* 


„ » " ° „ x 




0.4 


" 


B B 




0.2 




%o = 









+ " " + 


■ 


0.2 




. * , , 


" 


0.4 




+ + 


" 






cooling sieps 


+ 


0.6 




x jk * x 30 cooling steps 


X 






+ + 60 cooling steps 


X 


0.8 




+. x x + 90 cooling steps 

+ + 120 cooling steps 

4- x -i- 150 cooling steps 


H 



(x,20,20) 

FIG. 7. Profile functions of the Polyakov loops through the 
lattice (varying x at y = 20 and z = 20) during cooling. 



The Polyakov loops around the central cube form a 
hedgehog. They reflect the color directions of the 
link variables. A parallel-transport to the above color 
frame leads to parallel Polyakov lines in color direction 
( — 1,-1,-1). The profile function of the Polyakov loops 
through the lattice during cooling is shown in Fig. [7] 
and the course of the central Polyakov line is plotted 
in Fig. [5] We find two transitions, one smooth between 
cooling steps 70 and 90 and a jump after 130 cooling 
steps, where the monopole finally vanishes. The color 
directions of link-, plaquette- and Polyakov matrices of 
the central cube of the negative spherical vortex after 120 
cooling steps are depicted in Fig. 9. 
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FIG. 9: The arrows indicate the rotational vectors of link-, plaquette- and Polyakov-matrices around the central cube of the 
negative spherical vortex after 120 cooling steps. Link-vectors rotate around the center, color fluxes through plaquettes are 
aligned parallel in the ( — 1, — 1, — Indirection and Polyakov loops form a hedgehog. 



VII. CONCLUSIONS 



We reported on violations of the lattice index theorem 
for smooth, "admissible" gauge configurations of classi- 
cal, spherical center vortices, for both, overlap and asq- 
tad staggered fermions in the fundamental and adjoint 
representations. Numerically, the discrepancy equals the 
winding number of the spheres when they are regarded 
as maps S 3 -> SU{2). 

During cooling a spherical vortex on a 40 3 x 2 lattice, 
the index in the adjoint representation indicates evidence 
for an object with fractional topological charge, which 
is identified as a Dirac monopole with the well known 



singularity in its center and a gauge field fading away 
at large distances. Therefore even for periodic boundary 
conditions it does not need an antimonopolc. 
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